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Abstract 

The production of the isoscalar J 71 " = 3 + didelta dibaryon d* by proton 

inelastic scattering from deuteron targets is described in a double-scattering 

Glauber approximation. Each scattering changes a target nucleon into a A 

with the help of the isovector tensor force transmitted by tt and p mesons. 

The differential cross section constructed from empirical Love-Franey nucleon- 

nucleon i-matrices and a simple model of d* shows a maximum of some 10 

/xb/sr at 70° (cm.) for 500 MeV protons. The partial width of the decay 

d* — ► NN caused by the exchanges of the same mesons is found for this 

simple model of d* to be about 9 MeV if the d* mass is 2100 MeV. The 

implications of these results are discussed. 
PACS numbers: 14.20.Pt, 25.40.Ep, 23.50,+z 
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I. INTRODUCTION 



The successes of quantum chro mo dynamics (QCD) between quarks as the fundamental 
theory of strong interactions have led people to expect new hadronic states often dominated 
by exotic Fock-space components Jl|-0. In nuclear physics, one is particularly interested in 
new nonstrange dibaryons (with baryon number A = 2) of unusually low mass and narrow 
width that might betray their underlying quark structures. No such dibaryon has been 
unambiguously identified experimentally despite years of search 

There are two promising dibaryon candidates, one with unusually high mass and one 
with unusually low mass. The high-mass candidate is an isospin 1 structure of width not 
exceeding 80 MeV first seen experimentally in the helicity difference A<jl{pp) of the total 
pp cross section at an energy that corresponds to a dibaryon mass of 2735 MeV 0. It has 
been seen more recently at the same mass in the pp spin correlation parameter A 00nn J7|. 

This structure has been interpreted by Lomon and collaborators || el. S cl SIX - quark "small- 
bag" state with the nucleon-nucleon (NN) quantum number of 1 S'o. The interpretation 
uses the R-matrix formalism to determine if the matching of an internal quark description 
based on QCD to an external nucleon description based on meson-exchange dynamics at 
a boundary radius separating the two regions could be made in a way consistent with the 
empirical NN phase parameters in the neighborhood of the observed structure. (For an 
assumed quark model of internal wavefunctions, the resonance energy can be predicted 
by varying the matching radius tq until the external NN wavefunction vanishes at tq at 
precisely the same cm. energy as the internal bag-state energy.) The R-matrix analysis 
also yields a resonance width of about 50 MeV, in rough agreement with experiment. This 
interpretation will require confirmation by phase-shift analysis or by direct detection via 
resonance production in nuclear reactions. 

In the R-matrix analysis, the experimentally observed mass (2735 MeV) has been found 
|| to be consistent with the bag parameters used in the "Cloudy-Bag" model of ||, provided 
that the pion cloud is neglected. Furthermore, because of its coupling to the external NN 
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channel, the internal bag state is not in equilibrium, and therefore the resonance mass is 
higher than the equilibrium bag mass of 2680 MeV jn| . If the external pion cloud had been 



present, as in the actual Cloudy-Bag model 0, the equilibrium mass (at 2380 MeV) would 
have been much lower fTIfl . It thus appears that the role of the pion cloud external to the 
bag needs clarification. 

The R-matrix analyses have shown that NN phase parameters are also consistent with 
many other "high- mass" nonstrange dibaryons fl(| , such as those predicted by nonrelativistic 



potential models with pairwise color confinement It is obvious that the case for these 
high-mass dibaryons can be significantly strengthened if experimental effects are observed 
for at least another of these dibaryons. Experimental structures seen in Acr^Qop) at 2900 
MeV |L3| and in Aa^i^np) at 2630 MeV |L4|] could be candidates for a 3 Pi and a 3 Si dibaryon, 
respectively |L5 . 



In addition, suggestions have been made that these nonstrange dibaryons might appear 
at much lower masses instead. One with a proposed mass 2065 MeV and width T wNN = 0.5 
MeV might be responsible for a narrow structure in the energy dependence of the experi- 
mental excitation function [at 5° in the center of mass (cm.)] of the pionic double charge 
exchange (DCX) reaction nn(TT + , ir~)pp on nuclear neutrons at T n = 50 MeV |p.G(| . This 
explanation seems to be supported by observations of a narrow structure at 2060 MeV with 
a width < 15 MeV at ITEP |~7| and at CELSIUS | 18[| . The associated dibaryon, usually 



called d', has the proposed quantum numbers T = 0, J 77 = , making it inaccessible from 
NN channels and copnsequently narrow. However, an alternative explanation of the DCX 



phenomenon that requires no dibaryon has also been given [T9[ . 

The dibaryon interpretation finds theoretical support in bag models of dibaryon masses 
where this particular dibaryon appears at 2100 MeV |20] or 2000 MeV f2T| . The theoretical 
bag state involved is a P-wave excitation in the q 2 — q 4 separation with the cluster quantum 
numbers of (T, S)u = (0, 0) and (T, 5)3456 = (0, 1). However, in quark potential models with 



pairwise color confinement, the state appears much higher, at around 2700 MeV P2|. The 



mass can be reduced considerably with configuration mixing, but it seems difficult to reduce 



it to below 2400 MeV if one uses quark-quark (qq) dynamics deduced from single-baryon 
resonances 
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This paper is concerned with a second nonstrange isoscalar dibaryon called d* which 
has the quantum numbers J 71 " = 3 + , is accessible from np( 3 D^ — 3 G$) channels, and whose 
baryon-baryon component is made up of A 2 . We shall call this d* a didelta when we want 
to emphasize this A 2 component. 

The theoretical d* mass m* again covers a wide range: It is highest in the small-bag based 



R-matrix analysis which yields a value of 2840 MeV ||23|| , well above the A A threshold at 
2460 MeV. It is lowest in the Quark Derealization and Color Screening (QDCS) Model, 
where it appears at around 2100 MeV [HHH, just above the nNN threshold at 2020 MeV. 



Near the lower limit of this mass range, the irNN phase space is small so that the decay 
of d* is probably dominated there by the NN channel, where the nucleons fall apart in a 
relative D-state. However, the (ir) n NN widths could dominate as m* increases. 

A search for the d* dibaryon is interesting for the following reasons: In most quark 
models, the A — N mass difference comes from the color-magnetic term of the one-gluon 
interaction between pairs of quarks. The total pairwise color-magnetic operator has the 
same (repulsive) matrix element in d* as in two well separated A's if the former's orbital 
wavefunction is totally symmetric in the quark labels ||27|| . In the Massachusetts Institute 



of Technology (MIT) bag model f28 |, the d* mass then falls below the AA threshold by 120 
MeV because the spatial integral associated with the interaction is inversely proportional 
to the bag radius R, and this radius increases from A to d* by virtue of the increasing 
total kinetic energy |l|j24| . (The radius of the A-baryon bag of baryon number A is roughly 
proportional to A 1 / 3 in the MIT bag model.) 

The situation is different if quark confinement comes from a pairwise qq interaction that 
rises to infinity at infinite separation. The increasing size of an A-baryon containing 3A 
quarks causes its confinement energy to increase so much that the d* mass moves substan- 
tially above the A A threshold instead [|12||. This result has been confirmed by 



The QDCS model is able to reduce the d* mass substantially with the help of two 



additional assumptions: quark derealization (QD) and color screening (CS) |24|- [2"6|1 . QD 
takes advantage of the fact that the kinetic energy of a single quark state could be reduced if 
it is partly on the left side and partly on the right side of the system. For a Gaussian spatial 
wavef unction, a maximum reduction of the kinetic energy of about 19% appears when there 
is a 50/50 left/right separation with the two wavefunction centers separated by 2.3 oscillator 
lengths, like two peas in a pod. For a dibaryon built up of a product of six such delocalized 
quark wavefunctions, about 72% of the system is in the cluster configurations of the type 
q 2 — q 4 and q — q 5 where the reduction in the kinetic energy can be realized. If one reduces 
the total kinetic energy in the MIT bag by the resulting 14% (72% of 19%), one gets a d* 
mass of about 2090 MeV instead of the usual 2340 MeV, assuming that the interaction and 
confinement energies retain their spherical forms. 

However, in color confinement models or even in string models, this QD reduction of 
the kinetic energy alone cannot overcome the strong increase in pairwise color confinement 
energy with increasing qq separations and with increasing baryon number A. Color screening 
now comes into play by assuming that the rising repulsion of the confinement potential at 
large distances [I2j does not exceed a finite upper bound J23|,[^J. Now finally QD could 



become energetically favorable in /L-baryons. (The model of |2B| achieves the same result 
by measuring the confinement energy from the nearer baryon center, thus avoiding large qq 
separations.) 

The QD phenomenon could presumably be restricted to the interior of a bag. However, 
if realized between well-separated baryons |[25|| , the idea has far-reaching implications. With 



quark derealization taking place at all densities, the transition to the quark-gluon plasma 
will be at best a second-order phase transition. This seems to imply that if such an extreme 
picture is correct, attempts to search for quark-gluon plasmas might be doomed to failure, 
given the difficulty of detecting unambigous signals from even a first-order phase transition. 
However, a counter-argument is provided by the recent observation of unusually strong 
absorption of J/ip mesons in Pb-Pb collisions |3t]]. One interpretation is that this is a 



signal for the color deconfinement phase transition [p2]| . Consequently, it would be interesting 



to look for experimental indications for or against the QDCS model. This cannot be done 
by studying nuclear forces, which can already be understood in terms of meson exchanges. 
In contrast, the observation of a low d* mass could be taken as a signal for QD in A-baryons. 

In this connection, it is worth noting that a recent R-matrix analysis of available NN 
phase shifts below the dibaryon mass of 2240 MeV finds no sign of a d* resonance in the 
np( 3 D^ — 3 G3) channels with a width greater than 1 MeV [15|. However, the analysis does 



not exclude a narrower d* at a mass in between the energies of known phase shifts. 

In any case, progress in our understanding of dibaryons will require new experimental 
inputs. In particular, any new information on whether the d* mass might be high or low is 
likely to have important implications on the dynamics of quark confinement in baryons. 

The dibaryon d* could be produced by the inelastic scattering of projectiles from nuclear 
targets. The understanding of past failures to find it and the justification for future 
searches in such reactions would require some theoretical input concerning its production 
and decay properties. The main purpose of this paper is to explore how its production 
cross section in pd inelastic scattering and its partial decay width into two nucleons could 
be calculated using standard techniques in nuclear reactions and a simple model of the d* 
as a didelta object. Some of the issues which must be resolved before realistic results can 
be obtained are briefly discussed. 

II. THE PD -> PD* PRODUCTION CROSS SECTION 

In our treatment of the inelastic production, the deuteron is described by an S-state wave 



function made up of a sum of three Gaussians fitted to the Bonn C S-state wave function [33 
renormalized back to 100%. The d* is taken to be a pure A 2 Gaussian wave function with 



an average AA separation of 2r* = 1.4 fm [p4]-p6|. The quark wave function in each baryon 



is assumed to be the same in both N and A [Q. The quarks in each baryon are localized 



to the left or right A in d*, with no left-right antisymmetrization or derealization [23j] yet 
included. Later in the paper, we shall estimate qualitatively the effects of derealization in 
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the calculated quantities. 

The excitation of d(T = 0, 1 + ) to d*(T = 0, 3 + ) requires an isoscalar transfer of 2 units 
of angular momentum. In the usual models of nuclear forces containing the exchanges of 
only pseudoscalar, scalar and vector mesons, d* can only be reached in the lowest order by 
a spin-isospin flip in each of the two target nucleons. Each spin-isospin flip is caused by 
the exchange of an isovector meson such as ir and p between the projectile and the target 
nucleon. It will turn out that the unpolarized cross section has important contributions 
from that part of the NN t-matrix proportional to the operator (<7"i.q)(o"2.q)(Ti.T2). 

The spin- averaged differential cross section in the cm. frame has the structure 

da fi _ P*P* f da h 



dn* it dt 



P *p}(|A| 2 )spin, (1) 



where p* is the proton momentum in the reaction cm. frame in the initial or final reaction 
state a. The invariant inelastic amplitude 



A(q)«(^)/S fi (q,q') 



/i(qi) 72(02) 



d 2 q' (2) 



is approximated by the Glauber double-scattering contribution between quarks in the 
cm. frame of the reaction. Each scattering changes a target nucleon into a A. The integral 
involves two momentum transfers 

qi = |q + q' , q2 = |q - q' , (3) 

and an inelastic formfactor Sfi(q, q'). The actual expression is considerably more complicated 
than the symbolic structure shown in Eq. (|2]), as we shall discuss below. 

The outside factor in Eq. (0) contains the Glauber pair factor iV = N^Np, where iVr = 9 
for the number of distinct quark pairs on the target side, with one quark from each of the 
two target baryons. The effective projectile pair number Np is made up of a contribution 
of 3 from three projectile quarks each interacting twice with the target, and a contribution 
from three pairs of projectile quarks of 

3<±(<r 1 .ir 2 )(T 1 .T S ,)) = f. (4) 
7 



The formfactors from the projectile quarks are different for these two cases, being both 
different from those for two single scatterings from two separate baryons. The difference 
is described by an extra projectile formfactor Sp(q, q') above and beyond those for two 
separate elastic NN scatterings. Its functional form will be given later. 

There are quark-quark operators hidden in the quark-quark scattering amplitudes in 
Eq. (0). They are treated by first expressing NN scattering amplitudes in quarks coordinates 
pEfl . After the calculation of various operator matrix elements, the scattering amplitudes 
are re-constructed back to iViV form. 

The momentum k* is the NN relative momentum in the NN cm. frame. For elastic 
scattering, or in the high-energy limit where the inelasticity is negligibly small, it is sufficient 
to use the elastic scattering value k e \ = (3/4)p*. We recognize that in inelastic scattering 
at lower energies, the effect of the smaller projectile momentum p*f in the final state should 
be taken into consideration, in order to describe more accurately the energy dependence of 
both kinematics and dynamics. One possibility is to use the geometrical mean momentum 
(3/4)^Jp*p*j:, but this cannot be correct because it gives the wrong behavior at threshold. 
Of course, the Glauber multiple-diffraction formalism is a high-energy approximation that 
should not be used too close to a reaction threshold, but it is conceptually also important to 
ensure that the invariant amplitude A does not have a spurious singularity at threshold. The 
threshold behavior is not just a question of kinematics, because the elementary amplitude 
/i(q) is also involved. 

Now the inelastic production of d* involves at least two collisions. At each collision, 
the energy transfer can have a range of values. This suggests that the correct description 
might require an average over some distribution. We shall use the simplest realization of 
this concept, namely the assumption that the NN dynamics in Eq. (0) should be that for 
the arithmetical average of the proton energies in the initial and final states: 

^av = T lab - I Am*, (5) 

where Am* = m* — is the mass transfer in the inelastic scattering. The momentum k* is 
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the relative momentum in the NN cm. frame when a proton projectile of kinetic energy T av 
is incident on a nucleon target. In this prescription, the Glauber double-scattering integral 
for the invariant amplitude does not show a spurious singularity at threshold. 

The dynamical factors (fi/k*) inside the Glauber double-scattering integral is approxi- 
mately frame invariant, because in the usual parametrization using the optical theorem, it 
depends primarily on the total NN cross section <tnn if nucleons were the elementary objects 
of the Glauber method. 

The reduction of NN amplitudes involving the operator (<Ti.q)(cr 2 .q)(Ti.T2) to quark- 
quark amplitudes is made in a number of steps: 



a) We first use the well-known relation |35| g mqq = (3/5)g m NN) where m is 7r or p. 



(b) The qq operators are simplified by using the identities 

(<T.qi)(o\q 2 ) = qi.q 2 + i<r.(qi x q 2 ) , (6) 



(Oi.qi)(0>q2) = -(oi.o^Xqi.qa) + |(o-< x Oy).(qi x q 2 ) + (o- x ^^.(qi x q 2 ) (2) . (7) 

Only the tensor-force terms involving (cr 6 x (Tg) ( - 2 \Te-Tg) for two target quarks 6 and 9 
contribute to d* production. (Quarks 4-6 are in the first target baryon, while 7-9 are in 
the second.) On the projectile side, we keep only the term that is spin- independent. A 
projectile spin-orbit term can be shown not to contribute to the spin-averaged production 
cross section under rather general circumstances. 

The identity fl7|) can also be used to handle central spin-isospin dependent interactions 
by replacing one or more of the momentum transfers q m from the tensor interactions by 
projectile quark spin operators <r m of the central interactions. In this way, one can show 
that when both interactions are central, the production of d* induces 2 untis of spin transfer 
at the projectile, which cannot then remain a proton. However, the production pd — > pd* 
is possible when one interaction is central and the other is tensor, because in this case the 
projectile suffers only 1 unit of spin transfer. This mixed tensor-central contribution will 
not be included in the present exploratory study. 
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(c) In the result reported here, we have neglected all quark-exchange terms between 
target baryons. This permits the required operator matrix elements to be calculated in 
terms of the reduced matrix element 

((N 2 )S'T' = 10 || (<x 6 x <t 9 )( 2 )(t 6 .t 9 ) I) (A 2 )30) = fy/1. (8) 

Results with fully antisymmetrized quark wave functions for the target will be reported 
elsewhere when completed. 

On the other hand, the use of NN t-matrices ensures that all quark exchanges between 
projectile and target nucleons are automatically included. 

(d) The spin-orbit terms in the NN amplitude are neglected, as we are not concerned 
here with polarization phenomena. 

The spin-averaged differential production cross section then takes the form 



dtt* V15 



where 



^(q) = (|) / dVMq.qOBCq * q) - (q' * q')fl, (io) 
Mq,q') = (b 2 ~ i 12 ) (£) tU^)tU^)SMM)s f M), (12) 

and m is the nucleon mass. 

We shall use the empirical ./V./V isovector-tensor (IVT) t-matrix t- lvt constructed from AW 
phase shifts by Franey and Love (LF) P5| . The one appearing here is actually three times 
the tensor-force function tabulated by LF, i.e. the iViV scattering amplitude is here defined 



as 



/nn(q) = - (7-) *nn(q) 



47T/ 

m 



An 



*ivt(9)(o , i.q)(o-2-q)(Ti.T 2 ) + (13) 
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where the tensor operator is 1/3 of that contained in the usual tensor- force operator used 
in nuclear physics. 

In the Born approximation, £nn(<i) simplifies to the Fourier transform V(q) of the NN 
potential in the notation of Ref. ]33|]. In fact, we shall also use the Full Bonn potential |33 



in the Born approximation to study the separate contributions of it and p exchanges and 
the dependence of the production cross section on the projectile energy. 
The additional projectile formfactor needed in Eq. fll2l ) is 



Sp(q,q) = exp 



q 2 + f_ 



12b 2 3b 2 



exp 



q 
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, for single quarks; 
for quark pairs. 



(14) 



246 2 66 2 _ 

Each of these goes with its own effective projectile pair number, a minor complication that 
will not be shown in the formulas. 

On the target side, there are two separate scatterings with two baryons. When the 
inelastic scattering leaves the size of the target baryons unchanged, as we have assumed in 
this paper, the baryon formfactors themselves could be reabsorbed back into NN scattering 
amplitudes as they are reconstructed from the quark-quark amplitudes. We are then left 
with just the baryon-baryon (or wavefunction) inelastic formfactor 



Sfi(q) 



2(3(3* 



3/2 
2 / eX P 



q 



2((3 2 + (3* 2 ) 



(15) 



\P 2 + P* 

if the deuteron wavefunction is also a single Gaussian, with the inverse range (3, while 
(3* = J (3/8)/r* describes the d* wavefunction of r.m.s radius r*. (This wavefunction radius 
does not include the contributions from the baryon formfactors.) In reality, the Bonn C 
deuteron S-state wavefunction is used after being fitted to the nonorthogonal three-term 
oscillator form 



V'Bonnc(p) ~ ^ C ^G°)' 



(16) 



i=l 



where ipi{p) is a normalized S-state oscillator wave function. The resulting range parameters 
obtained by minimizing the percentage m.s. deviation are 
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7 = (7i, 72, 7s) = 0.04467 (1, 5.04, 21.5) fm 



-2 



(17) 



where j { = 2(3?. The expansion coefficients, renormalized from the fitted value of 94.34% 
back to 100%, are 



The momentum transfer q of the reaction is of course calculated with the correct (rel- 
ativistic) inelastic kinematics. However, the other momentum transfers of the Glauber 
approximation are treated in the high-energy limit where the inelasticity is negligibly small. 
As is known, the z-axis of the Glauber formula is usually not constant in space but cho- 
sen instead along the bisector of the initial- and final-state momenta. This means that the 
momentum transfer is on the equatorial, or xy, plane. In fact, all momenta in the Glauber 
double-scattering formula lie on this equatorial plane. In extending the formula to inelastic 
scattering, we have kept all momenta on the equatorial plane even when the energies in- 
volved are not very high. This means that two of the terms in the sum in Eq. (0), namely 
for jj, = ±1, are zero. 

III. RESULTS FOR THE PRODUCTION CROSS SECTION 

The results shown in this section are all calculated in that angle-averaged approximation 
in which cos#' everywhere inside the integrand in Eq. (|2|) is taken to be l/y/2, 9' being the 
angle between q and q'. The effect of a full angle integration will be discussed near the 
end of this section. For the sake of completeness, we shall give full angular distributions 
even though the Glauber multiple-diffraction approximation is known to be reliable only for 
small angles. 

Fig. 1 gives the differential production cross sections at the TRIUMF proton energy 
of 516 MeV. The 1985 Love-Franey (LF) t-matrix [|36] at the lab energy of 515 MeV (on a 
nucleon target) is used for d* masses of m* = 2050, 2100 and 2150 MeV, a range of particular 
interest for the QDCS model mentioned in the Introduction. Calculated with the same NN 



c = (ci,c 2 ,c 3 ) = (0.31491,0.49716,0.36926). 



(18) 
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t matrix, these cross sections involve exactly the same invariant production amplitude. They 
differ only in the outside kinematical factors p*p*f in Eq.([L]) and in the fact that the angular 
distribution covers a narrower range in momentum transfer as m* increases. For this reason, 
it is often necessary to show results for only the smallest m* = 2050 MeV because it has the 
largest range of momentum transfers. 

The cross sections shown are about 7 /xb/sr at the second maximum at about 70° cm.. 
They decrease with increasing d* mass, and should vanish at threshold. The structure 
at small angles, seen very clearly in the m* = 2050 MeV result, comes from the /i = 
component of the production tensor operator. Its contribution to the differential cross 
section of Eq. (|9|) for m* = 2050 MeV is shown separately as a short dashed curve in fig. 1. 
The remaining contributions are from the \x = ±2 components. Being proportional to g 4 for 
small q 2 , they are responsible for the second maximum. The strong fi dependence means that 
the angular distributions can be expected to be very different when the colliding particles 
are polarized. 

The production mechanism, requiring at least a double scattering, is particularly sensitive 
to the dynamics of the effective qq tensor force, here derived from the empirical NN tensor 
force. It is therefore of interest to show, in Fig 2, how strongly the production cross section 
for m* = 2050 MeV at the TRIUMF energy increases as the effective nucleon lab energy of 
the input NN t-matrix is decreased. 

Since the production cross section as calculated depends only on the isovector tensor 
part of the NN i-matrix, the effect seen is a direct reflection of the latter's rapid increase 
with decreasing nucleon energy. This behavior is well understood in the theory of nuclear 
forces: The shorter-range part of the tensor force due to the exchange of p mesons is opposite 
in sign to the longer-range part from 7r exchange. As the scattering energy decreases, the 
strong and long-range 7r-exchange contribution becomes rapidly more dominant. We shall 
be able to show separately the results calculated from each of these two parts of the NN 
tensor force when we use the Full Bonn potential. 

The inelasticity involved in d* production for the m* range examined here is a very 
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appreciable fraction of the energy available in the beam. Given the strong dependence of 
the production cross section on the effective nucleon energy of the NN isovector tensor 
i-matrix, it is necessary to choose this energy carefully. For reasons discussed in the last 
section, we shall use the average lab energy T av shown in Eq. (|5[) which is the arithmetical 
average over the initial and final states. 

The resulting cross sections for different d* masses are shown in Fig. 3. These are 
obtained by interpolating the results calculated for the three lab energies T = 515, 425 and 
325 MeV tabulated by LF. We see that interpolated results have increased by a factor of 
about 2.0 over the value calculated at the incident energy for in* = 2100 MeV. The final cross 
sections at the second maximum are about 13 /xb/sr. This is very large, but it is perhaps 
not totally unexpected because the production involves the same amplitudes responsible for 
the resonance production of A from nucleon targets. 

The use of empirical t-matrices takes care of re-scattering or wave-distortion effects in 
the sense of an impulse approximation. To determine how important such effects are, we 
use the Full Bonn potential treated in the Born approximation. The results obtained with 
the NN relative momentum k* shown in Eq. (|2p calculated from T av are shown in Fig. 4 for 
different values of m* . 

We see that these production cross sections have the same angular behavior as those in 
Fig. 3 for the LF t-matrix, but their values are higher by a factor of about 2.8. The effect 
might seem huge, but since the calculated cross section is proportional to g 8 , where g is a 
meson- iViV coupling constant, the t-matrix reduction of the effective g 2 is by a very modest 
factor of 0.77. In other words, a strong sensitivity to NN dynamics is unavoidable in such 
a high-order production process. A quantitative calculation of the cross section might well 
require better t-matrices constructed from more recent iViV phase shifts. 

The use of the Bonn potential allows the contributions of the n- and p-exchange poten- 
tials to be separated. The cross sections for each contribution alone are also shown in fig. 4. 
Note that their total contribution to the cross section is not the sum of their separate con- 
tributions because amplitudes interfere and the production mechanism is double scattering. 
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The 7r-only result is enormous, but it is effectively controlled by the much weaker p ex- 
change contribution. Since this cancellation is not a well-determined part of NN dynamics, 
Fig. 4 also shows the importance of using the best NN input in the calculation. 

The Born approximation should improve in accuracy with increasing projectile energy. 
Fig. 5 shows how the calculated cross section for m* = 2050 MeV decreases with increasing 
energy, while the second maximum moves forward in the angular distribution, but not in 
the momentum transfer. Both features are consistent with the expected energy dependence 
of cross sections. Of course, the accuracy of the Bonn potential at these higher energies is 
somewhat uncertain. In spite of this reservation, Fig. 5 does suggest that the production 
cross section will not decrease sharply with increasing energy. 

The production cross section also depends on the d* wavefunction size r* (half of the 
average AA separation). Fig. 6 gives the results calculated for r* = 0.5, 0.7, 0.9 fm, covering 
a realistic range of possible d* sizes. Other parameters used are m* = 2050 MeV and an 
effective nucleon energy of 425 MeV for the LF t-matrix, very close to the recommended 
average value of 429 MeV. We see that the effect is fairly strong especially for smaller 
d*'s presumably because the momentum transfers involved can then be more different. We 
therefore conclude that the calculated cross section can be quite sensitive to the short-range 
components of the wave functions of both d and d*. The short-range components neglected 
in the present calculation include the deuteron D-state, exotic admixtures, and the effects 
of quark exchange and derealization. 

Finally, we address the question of the accuracy of the angle- averaged approximation 
used this section. A selected number of angle-integrated cross sections for the TRIUMF 
energy using the Bonn potential with m* = 2050 MeV have been calculated. The angle- 
integrated result at the CM angle 6* = 70(0,180)° is 40.31 (31.06, 6.247) /ib/sr when the 
angle-averaged approximation gives 38.54 (32.06, 6.252) /ib/sr instead. Thus the angle- 
averaged approximation seems to have adequate accuracy. 

The present calculations have given only a crude picture of the inelastic production of 
d*. Their many limitations will be discussed in the last section. 
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IV. THE D* -> NN DECAY WIDTH 



The dibaryon d* cannot decay into two nucleons if its constituents do not interact. The 
simplest interaction that can do it is a two-quark interaction containing the spin (<Tj x cr,-)^, 
to bridge the gap between initial and final intrinsic spins. If d* is a didelta, even when these 
deltas are not pointlike, the perturbing operator must also contain the isospin operator Tj.Tj 
in order to change a A into an N at each quark vertex. It is clear that isospin-independent 
interactions cannot do it. Thus all direct gluon exchanges do not contribute, no matter how 
strong they are. 

Among the simplest two-quark interactions that can do it are the one-meson-exchange 
potentials carried by tt and p. We shall show below that their contributions are quite large. 

It can be argued that these mesons might not have the presence inside d* that they 
have in the outside meson cloud, because much of the dibaryon interior might be in a 
different vacuum state, the so-called perturbative QCD vacuum, where mesons loose their 
individuality. Since QCD is flavor- independent, the required Tj.T, operator could only come 
from the Heisenberg isospin exchange operator 

V T = |(1 + T i .T i ) (19) 

arising from the Pauli exchange of two quarks. We now show that the resulting potential 
has a functional form similar to that for one-meson exchange. In this preliminary study, we 
shall not consider any three-quark interactions, including those where two interacting quarks 
involve a noninteracting quark via Pauli exchange, because they have more complicated 
structures. 

The exchange two-quark interaction generated by one-gluon exchange (OGE) between 
quarks i, j can be written in the familiar form p7tp8| 

V xm = -V l3 V m = -V X V T V°V X V^ (20) 

where 
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V X = \ (1 + VA,) (21) 

is the color exchange operator. The space exchange operator V x interchanges the spatial 
labels in the qq final state. For an exchanged gluon of effective mass /i, taken here to be 300 
MeV [ftm , the direct qq interaction from OGE has the standard one-boson-exchange form 

V M = 7tt^( a - a ^) \-^M 2 + (^-q)K--q) + •••] • (22) 

The QCD coupling constant appearing in Vo is dependent on the gluon mass fi. In addition, 
it should be an effective coupling constant that might include certain higher-order effects. 

Of the terms shown, the second term is a tensor interaction which tends to be unobtrusive 
in baryon spectroscopy, just like tensor forces in nuclear spectroscopy. However, the first 
term is the color-magnetic interaction which can be related to the A — N mass difference 
Am = m A -m N = 293 MeV: 

< 23 > 

Here 

= < 24 > 

is the spatial two-quark matrix element in a baryon (B), taken for simplicity to be the same 
in both A and N. There are momentum-conserving 5-functions 5k for the two quarks i,j at 
each gluon-quark-quark vertex. 

Is this empirical qq interaction strong or weak? This question can be answered by 
anticipating that even though both V and the gluon propagator in V qq depend significantly 
on the gluon mass fj,, their effects tend to cancel when the same A — N mass difference is 
fitted. The interaction strength Vq turns out to be 0.11 GeV -2 when \x is taken to be the p 
meson mass. In contrast, the equivalent strength of the one-rho-exchange (ORE) potential 
in the Full Bonn potential is 
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MD^fH 1+ tf^ 54GeV ~- (25) 

where the factor (3/5) 2 comes from the reduction from NN to qq operators. This is about 
500 times stronger than the qq interaction strength. Although the pNN coupling constant 
is not well determined in the NN interaction, it is clear that the ORE potential appropriate 
to the baryon exterior is some two orders of magnitude stronger than the effective qq inter- 
action present in the perturbative vacuum of the baryon interior. This shows that the OGE 
contribution to the partial decay width is negligibly small compared to the meson-exchange 
contributions. 

It is nevertheless interesting to complete the derivation of the exchange qq interaction 
and to determine how the space-exchange operator V x further affects the final result. We 
therefore go on by noting that the term in V^ qq containing the d* — > NN decay operator 
can be isolated by using the expansion 

- V x V T V^\ l .X 3 )( ( T l . q )( ( T rq ) = -|(«r i .q)(<r i .q)(T i .T i ) + ... , (26) 

obtained with the help of the identities 

(A..A,) 2 = f - KA..A,) (27) 

and 

(cr i .cr i )(cr i .q)(cr j .q) = (cr i .q)(cr j .q) + (1 - (Ti.er^q 2 . (28) 
The final result is 



K qq (q) = ~V X ( - ) (<r i .q)(*r i .q)(T 4 .T i )- J -^ + ... , (29) 



showing only the term which can turn colorless A's into colorless iV's. This term too has a 
strength characterized by Vq. This is the only term that can contribute, in the lowest order, 
to the decay of the d* treated didelta. 

In more realistic models of d*, colored A's appear in the so-called "hidden-color" com- 



ponents |]40 |. These colored objects can decay into iViV by both direct and exchange OGE 



potentials via terms proportional to the color operator Aj.Aj. However, all such contribu- 
tions are necessarily based on the effective qq interaction V qq whose strength V is two orders 
of magnitude weaker than meson-exchange interactions. 

We are now in a position to calculate the decay width. Its spin-averaged value in first- 
order perturbation theory EIJ has the structure 



r(m*) = / d 2 n p ,\N m ((N 2 )p*\V\(A 2 )d*)\% in ^ (30) 

where p* is a nucleon momentum in the center-of-mass frame and fi*j (= m*/4) is the 
relativistic reduced mass in the final state calculated with dynamical masses or total energies. 
iV qq = 9 pairs of qq interaction V contribute to the decay. The momentum state |p*) has 
the normalization: 

After some algebra, the decay width can be written in the final form 



1 m = 



qq 9 V 2 



15 (27T) 2 7T 3 / 2 V 

where the model-dependent factor F mo d e i is 



16 1 ((3 * 7 \f (m\ 

— r model; l^J 



Fmx — (|) 2 Iivt( K i lj 1) 1) (33) 



for the meson-exchange (MX) model, and is 

gp,2 \ 3/4 



f -=UJUf^J w«,A.ft.o.) (34) 

for the exchange-gg (XQQ) model. Here 

I mo U^D u D 2 ,D 3 ) = e - D ^ 2/2 J e~ D2Q2/2 [j 2 (iD 3 KQ)(KQ) 3 }t modcl (P*Q)QdQ, (35) 

k = k*//3* is the dimensionless nucleon momentum in the cm. frame, and Q = q//3* is the 
dimensionless momentum transfer. The elastic baryon form factor at each quark vertex is 
already contained in the NN interaction t- lv t of Eq. flT3"|). Furthermore, the t-matrix contains 
iViV rescattering contributions to all orders of the NN interaction. 
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For the XQQ model, for which i xqq in this equation refers to the function Vo/(/i 2 + q 2 ) in 
Eq. (p9|), the exchange baryon formfactor appears explicitly in F xqq . The latter depends on 
the dimensionless parameters Di which are functions of the baryon size parameters (3* and 
b = 1/Vp, where r v = 0.6 fm is the proton radius. For r* = 0.7 fm, they have the numerical 
values of 

L>i = 0.214; D 2 = 0.398; D 3 = 0.153. (36) 

V. RESULTS FOR THE PARTIAL DECAY WIDTH 

The partial decay width for d* — > NN can now be calculated numerically for the empirical 
Love-Franey NN t-matrix, the Full Bonn potential treated in the Born approximation, and 
the OGE qq interaction with Pauli exchange. 

The empirical LF t-matrix used is evaluated at the equivalent nucleon lab energy 

Tfl* 2 

Ti a b = - 2m (37) 

2m 

of the final NN state. The resulting decay widths T, obtained by interpolating from the 
values calculated at the three nearest energies tabulated by LF, are shown in Fig. 7 as 
functions of the d* mass m* for different d* wavefunction radii r*. The dependence on these 
parameters is significant, but not strong, and is of a somewhat complex character. 

The origin of some of the complexities in the behavior of T can be seen in Fig. 8, which 
gives the corresponding results for r* = 0.7 fm covering a wider range of dibaryon mass 
calculated with the Full Bonn potential treated in the Born approximation. The result is 
typically 2-3 times larger than those for the LF t-matrix in the mass range (2050-2130 MeV) 
of interest in the QDCS model, in qualitative agreement with the corresponding behavior 
in the production cross section. Note however that the leading decay process studied here 
is first-order in the interaction, but the leading production process is second order. 

The "7r only" contribution (shown as a dash-dot curve) is again much larger, here by 
another factor of three in this mass range. Even the "p only" result (given by the short 
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dashed curve) is large, becoming in fact arger than even the 'V only" contribution after 
2690 MeV. This behavior is consistent with the fact that the decay involves interactions 
at shorter distances and higher energies than the production at TRIUMF energies. This 
means that the results for V is particularly sensitive to the cancellation between the n and 
p exchange interactions in our simple leading-order model. In fact, the vanishing of the 
calculated decay width due to the complete destructive interference between the it and p 
contributions can be seen in the figure at 2690 MeV. 

For the XQQ interaction appropriate to the perturbative vacuum of the baryon interior, 
we find a result of F = 77 eV at r* = 0.7 fm, m* = 2100 MeV and a gluon mass of 
/i = 300 MeV if V xqq interaction of eq. (|29"D is treated as a direct qq interaction. (This is 
done by dropping the space-exchange operator V x and adding baryon elastic formfactors to 
the expression for T.) Compared to the result of 22 MeV for "p only", this is a factor of 
3 x 10 5 smaller, in good agreement with the estimate based on the interaction strengths. 

On restoring the space-exchange operator V x , we find that the width drops drastically 
by another factor of 3 x 10 6 to only 2.1 x 10~ 5 eV. The result is very sensitive to r*, being 
250 times larger for the smaller d* radius of r* = 0.5 fm. One reason why the result is 
so small is that the orbital angular momentum in the relative baryon-baryon coordinate is 
L = 2, where the centrifugal potential inhibits the quark exchange. If we had used L = 
instead, the result would have been 2.3 x 10~ 3 eV, as compared to an L = width of 190 
eV when V x is purposely dropped. We conclude from this that in the absence of significant 
short-range components in the wave function, quark-exchange effects are negligible. 

Up to this point, we have done the calculation as if the interaction were 100% by meson 
exchanges via the meson cloud in the baryon exterior, or 100% by gluon exchange in the 
perturbative vacuum of the baryon interior. In reality, there is always some "external" 
contribution even when all the six quarks of d* are in the same cluster. Since the external 
contribution is so much greater than the internal contribution, some of it will always survive 
to give a nonnegligible decay width. In the absence of a detailed model describing the precise 
proportion between the two contributions, we shall discount our calculated meson-exchange 
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contribution by 50% as a guess of what a more realistic width should be. 

The interior correction for gluon exchange discussed here does not apply to the lowest- 
order production calculated in Sec. 3: The two mesons whose exchanges between projectile 
and target nucleons are responsible for the inelastic production of d* are both external virtual 
mesons. 

Is the meson-exchange interactions used here too strong? This can be answered partially 
by using the same method to calculate the decay width of that archtypical baryon decay 
A — > ttN. Using the ttNN coupling constant of /^nn/4tt = 0.078 given by the Full Bonn 
potential and reducing to a irqq coupling constant, we find a result of only 70 MeV, much less 
than the experimental value of 120 MeV. This result is in agreement with past calculations 
of this width fl35| . One might be tempted to increase the calculated result by a factor of 
(120/70) 2 ~ 3, but this is not advisable because a real 7r on the energy shell is emitted in the 



A decay, whereas the d* — > NN decay involves virtual mesons off the energy shell j|2) . The 
off-shell coupling constants that appear are more appropriately determined from nuclear 
forces, at least in principle. If we had been calculating the d* — > ttNN decay by the same 
method, we would be justified to increase the calculated result by a factor 120/70 ~ 1.7 for 
that vertex emitting the real pion. 



VI. DEPENDENCE ON QUARK MODELS 

The calculated d* production cross section and decay width can be expected to depend 
sensitively on the quark wavefunction of the d*, perhaps even more so than its theoretical 
mass. We shall consider qualitatively some of the issues involved. 

One of these issues is the possibility of quark derealization, which refers to the idea that 
under certain circumstances a quark may find it energetically favorable to be partly on the 
left side and partly on the right side of a dibaryon. The QDCS model actually describes 
each quark in the d* wavefunction as 50/50 left or right. The six-quark d* wavefunction 
then has the structure 
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(L + R) 6 = L 6 + QL 5 R + 15L A R 2 + 20L 3 R 3 + 15L 2 R 4 + QLR 5 + R 6 

-> 2(L 6 orP 6 ) + 12(L 5 RorLR 5 ) + 30(L 4 P 2 or L 2 P 4 ) + 20(L 3 R 3 ). (38) 

After the projection of relative S-states between the clusters and a correction for the position 
of the center of mass, the wavefunction simplifies to the form shown after the right arrow. 
Its components fall roughtly into two groups: There is a group of normal ( "n" ) clusters of 
q 3m configurations made up of L 3 R 3 , L 6 and i? 6 with no delocalized quark. They have the 
probability of 

P n = (2 2 + 20 2 )/1448 « 0.28. (39) 

The remaining group of components L 5 R,LR 5 ,L A R 2 , and L 2 R A has one delocalized quark 
("dq") away from a normal q 3m configuration, and the remaining probability of Pd q = 
1 - P n ~ 0.72. 

For the normal group, the projection of S-states makes the wavefunction spherical sym- 
metric in the relative baryon-baryon coordinate, and very similar to the Gaussian wavefunc- 
tion of our didelta model. In fact, the maximum overlap between the two wavefunctions is 
close to 100%. There is to be sure some depression of the two-center relative wavefunction 
near the origin of the relative coordinate, but the effect is quite unimportant in the wave- 
function overlap. The behavior of the short-distance wavefunction is probably much more 
important in the production and decay processes considered in this paper, but it is likely 
that the short- distance wavefunction is not very good in both models. Furthermore, the q 6 
component is entirely absent in the didelta model and is probably too weak in the QDCS 
model. 

Though subject to these additional uncertainties, our counting suggests that this normal 
group will contribute essentially the full amount, i.e. about 28% of that calculated in our 
model in both decay and production. 

For the abnormal components with one "wrong-way" quark, the contribution could be 
very different, especially if there is special coherence between the normal and abnormal 
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amplitudes. We are not in a position to estimate such coherent contributions because it 
would require a specific model. As far as its incoherent contribution is concerned, the worse 
that can happen is that it will vanish. This must be a rather extreme situation, because three 
of the nine pairs of interacting quarks involve the "wrong-way" quark, and the interaction 
could scatter it back to form a normal cluster structure. In the remaining six pairs, the 
"wrong-way" quark is a spectator, which requires a wavefunction overlap to get back to 
normal. There is thus a reduction in the calculated decay or production amplitude of the 
order of 1/2 or 1/e for the spectator contribution, more if the clusters are farther apart. We 
then end up with an estimate for the decay width or production cross section of the order 
of 7-20% from these abnormal components. 

Thus very crudely, we expect the derealization to reduce the calculated decay width or 
production cross section by a factor of 1/2-1/3. 

Another model-dependent issue is the contribution of hidden-color (HC) configurations. 
Our didelta model used without quark exchanges between the two baryon clusters contains 
no HC component. In contrast, most quark models of the d* contains 80% HC components 
where the first three quarks are in a color-octet state. These HC components are expected to 
contribute less, perhaps significantly less, than the baryon-baryon components. A calculation 
of their contributions using the method of this paper is now underway. For the time being, 
we shall allow for some contributions from the HC components by reducing the baryon- 
baryon results by a factor of 1/2. (The reduction factor is 1/5 when the HC components 
contribute nothing.) 

The final educated guesses for d* decay and production for the quark-delocalization model 
used with the Love-Franey NN t-matrix at m* — 2.1 GeV and r* = 0.7 fm are as follows: 
The decay width is decreased from 9 MeV to about 1 MeV when an "interior" correction 
of 1/2 and an octet-octet reduction factor of 1/2 are also included. The production cross 
section at TRIUMF energy is reduced from 13 to about 2 //b/sr at the second maximum. 
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VII. DISCUSSIONS AND CONCLUSIONS 



The leading-order processes studied here suggest that the inelastic production cross sec- 
tion of d* might be in the /xb/sr range, while its decay width into two nucleons might be 
in MeV's. Love-Franey empirical NN t-matrices are used to include all NN re-scattering 
effects to all orders in the decay and in an impulse approximation for the production. Other 
aspects of the reported calculations are not sufficiently realistic because of approximations 
made in the hadron wavefunctions and in the treatment of the reaction mechanisms. It is 
worthwhile to enumerate the most serious of these problems. 

The Glauber multiple-diffraction model used in the calculation of the production cross 
section might be quite good for elastic scattering at small angles. Its validity for large 
inelasticities and at large angles is unknown. It is necessary to correct for effects neglected 
by the Glauber model, especially at large angles ||43|| . However, there is probably no point 
in doing this unless one can also include higher-order production processes. These are the 
usual difficulties connected with the calculation of strong-interaction cross sections, and as 
usual, we see no simple solution. 

The Gaussian model of d* used here is very crude. Depending on the model, we need to 
add derealization and short- distance refinements. The S-state wave function used for the 
deuteron target is also inadequate. It is not difficult to put in the D-state NN component. 
Other short-range components such as hidden-color configurations can readily be treated 
too. The difficulty lies instead in the lack of knowledge on how strong these components 
are. Studies of the effects of short-distance wavefunctions require sustained efforts. 

In the calculation of the partial d* — > NN decay width, it is necessary to account more 
carefully for the role played by external meson exchanges versus internal gluon exchanges. 
It is obvious that this too cannot be done on a quantitative basis without a more realistic 
model of the d* wavefunction. It will be necessary to include higher-order processes not yet 
included by using NN t-matirces. 

In addition, the decay width could be dominated by the (Tc) n NN branches for sufficiently 
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large d* mass. It is necessary to understand these partial decay widths at least qualitatively. 

The theoretical picture concerning the dibaryon d* at the present time seems to be as 
follows: Its calculated mass has been in the range 2050 |24|]-2840 || MeV. Its inelastic 
production cross section could be significant, i.e. in the /xb/sr range. Its partial decay width 
into two nucleons is probably in MeV's for the low-mass candidate. Because its calculated 
mass is so sensitive to certain assumptions concerning quark dynamics in hadrons, any 
positive or negative experimental information on its presence in a certain mass range has 
interesting implications. 

What is the experimental situation concerning isoscalar dibaryons? A dibaryon search 
was made at Saturne by measuring the spectra for missing masses between 1.9 and 2.35 
GeV using the dd — > dX reaction for 2.29, 2.00 and 1.65 GeV deuteron beams. An upper 
limit of 30 nb/GeV 2 was found for the invariant production cross section of a dibaryon if 
its width is less than several tens of MeV [|J. This result is for the missing mass of 2.16 
GeV and for a 2.00 GeV deuteron beam with deuterons detected at 27° (lab), or 69.3° 
(cm.). It corresponds to a cm. differential production cross section at this angle of only 
15 nb/sr in the dd reaction. It is not easy to extract a pd bound from this result partly 
because of the presence of an elastic formfactor for the intact deuteron |44| , which causes a 
large reduction in the dd production amplitudes relative to the pd amplitude for production 
from single nucleons in the intact deuteron. An additional complication is that for the dd 
reaction, the d* can also be produced by another double-scattering process that involves 
both nucleons of the intact d. Its contribution can be expected to be similar in structure, 
but probably reduced in value, when compared to that of production from a single nucleon 
in the pd reaction. This process must also be included in the interpretation of dd cross 
sections. This means that any extraction of a pd bound from the dd bound will depend on 
a model-dependent theoretical analysis, and cannot be a pure experimental bound. 

To my best knowledge, the only direct experimental upper bound for resonance pro- 
duction in the pd reaction is an unpublished LAMPF experiment based on the d(p,p)X + 
reaction at T p = 798 MeV and 15.1° (lab). The results are in the range of 1-4 /ib/sr/MeV 
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H dependent on the missing mass in the missing-mass range of 1865-2200 MeV. They are of 
the same order as the very rough theoretical estimates of d* production given in this paper. 
Hence no definite conclusion can be drawn from a comparison between them. 

It thus appears that the present theoretical picture is still very unrealistic and incomplete. 
Much additional work is needed, especially on the partial decay widths in pionic channels 
for which there is at present very little quantitative information. However, the question of 
dibaryons is ultimately an experimental question. A new dibaryon search with a sensitivity 
much greater than the known LAMPF bound will be needed to advance our understanding 
of the situation. 
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FIGURES 

FIG. 1. Center-of-mass differential cross section for pd — > pd* at the proton lab energy of 516 
MeV for different d* masses to* using the 1985 Love-Franey t-matrix at 515 MeV. The to* = 2050 
MeV result from only the fi = term of Eq. (9) is also shown. 

FIG. 2. Differential cross sections for pd — ► pd* at 516 MeV for different Love-Franey t-matrix 
energies using to* = 2050 MeV. 

FIG. 3. Differential cross sections for pd — > pd* at 516 MeV for different to* masses using the 
Love-Franey t-matrix at the energy averaged over the initial and final states. 

FIG. 4. Differential cross sections for pd — > pcT at 516 MeV for different to* in the notation of 
Fig. 3 using the Full Bonn potential in the Born approximation with the NN relative momentum 
k* calculated at an average energy T av . The to* = 2050 MeV results for 7r-exchange only and for 
p exchange only are also given. 

FIG. 5. Differential cross sections for pd — > pd* for different projectile energies using to* = 2050 
MeV, the Full Bonn potential in the Born approximation, and an average nucleon energy T av . 

FIG. 6. Differential cross sections for pd — > pd* at 516 MeV for different d* wavefunction radii 
r* for to* = 2050 MeV using the Love-Franey t-matrix at 425 MeV. 

FIG. 7. Decay width for d* — > NN as a function of the d* masss to* for different d* wavefunc- 
tion radii r* using the Love-Franey i-matrix at the final-state energy. 

FIG. 8. Decay width for d* — ► NN as a function of the d* masss to* for r* = 0.7 fm using 
the Full Born potential in the Born approximation. The result for the Love-Franey t-matrix is also 
shown as a solid curve. 
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